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1 Introduction 



At first sight quantum chromodynamics (QCD) appears to be a relatively simple 
theory with extensive symmetry properties and only few parameters. To the novice 
it must come as a surprise that the whole range of strong interaction phenomena 
can be described by such a beautiful mathematical structure. Of course we do not 
know for sure whether this is really the case, but the accumulated evidence strongly 
suggests this to be so. 

The masses of the light hadrons and many other properties of these particles 
are very precisely known from experiment. On the theoretical side, however, it 
is difficult to compute these quantities from first principles, because of the non- 
linearities of the basic equations and because the coupling constant is not small at 
low energies. Precision tests of QCD in this regime have consequently been rare 
and are essentially limited to those aspects of the theory that are determined by 
the (approximate) chiral and flavour symmetries. It is evidently of great importance 
to overcome this deficit and to verify that the same lagrangian which describes the 
interactions between quarks and gluons at high energies also explains the spectrum 
of light hadrons and their properties. 

1.1 Lattice QCD 

To be able to address this problem one first of all needs a formulation of the theory 
which is mathematically well-defined at the non-perturbative level. Such a frame- 
work is obtained by introducing a space-time lattice and discretizing the fields and 
the action. The lattice cuts off the high frequencies and makes the theory com- 
pletely finite. Eventually the continuum limit should be taken and one then en- 
counters the usual ultra-violet divergencies which can be removed through the usual 
parameter and operator renormalizations. Lattice QCD may thus be regarded as 
a non-perturbative regularization of the theory with a momentum cutoff inversely 
proportional to the lattice spacing. 

Today quantitative results in lattice QCD are almost exclusively obtained using 
numerical simulations. Such calculations proceed by choosing a finite lattice which 
is sufficiently small that the quark and gluon fields can be stored in the memory of 
a computer. Through a Monte Carlo algorithm one then generates a representative 
ensemble of fields for the Feynman path integral and extracts the physical quantities 
from ensemble averages. Apart from statistical errors this method yields exact 
results for the chosen lattice parameters and is hence suitable for non-perturbative 
studies of QCD. 

One should however not conclude that the theory can be solved through nu- 
merical simulations alone. Without physical insight and ingenious computational 
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strategies it would be quite impossible to obtain accurate results that can be com- 
pared with experiment. 

1.2 Continuum limit 

One of the most important and difficult issues in lattice QCD is the continuum limit. 
In practice the limit is taken by computing the quantities of interest for several values 
of the lattice spacing a and extrapolating the results to a = 0. An obvious problem 
is then that one cannot afford to perform numerical simulations at arbitrarily small 
lattice spacings. In hadron mass calculations, for example, current lattice spacings 
are usually not much smaller than 0.1 fm. This will remain so for quite some time, 
because the simulation programs slow down proportionally to a 5 (or even a larger 
power of a) if all other parameters are held fixed. 

In view of these remarks, it is not totally obvious that the continuum limit can be 
reached in a reliable manner from the presently accessible range of lattice spacings. 
A detailed theoretical understanding of the approach to the limit is certainly required 
for this and extensive numerical studies are then needed to confirm (or disprove) the 
expected behaviour. 

1.3 Non-perturbative renormalization 

Non-perturbative renormalization is another topic where one would not get very far 
without theoretical input. Traditionally renormalization is discussed in the frame- 
work of perturbation theory, but when one is dealing with the low-energy sector 
of QCD the renormalization constants should be computed non-perturbatively. In 
more physical terms, the question is how precisely the perturbative regime of the 
theory is connected with the hadrons and their interactions. Evidently this is a 
fundamental aspect of QCD and not particularly a lattice issue. 

Non-perturbative renormalization involves large scale differences and is hence not 
easily approached using numerical simulations, where the range of physical distances 
covered by any single lattice is rather limited. The problem has attracted a lot of 
attention recently and various ways to solve it have been proposed. 

1.4 Contents & aims of the lectures 

The continuum limit and non-perturbative renormalization are the main themes of 
the lectures. In both areas significant progress has been made in the past few years. 
Chiral symmetry plays an important role in some of these developments and the 
way in which this symmetry is realized on the lattice is hence also discussed. 

Another item which is given considerable weight are finite-size techniques. Such 
methods are probably unavoidable if one would like to overcome the limitations 
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imposed by the currently accessible lattices. Some of them have already proved to 
be very powerful and the aim of the lectures partly is to illustrate this and to provide 
a better understanding of the underlying physical picture. 

The lectures are self-contained, but some familiarity with the basic concepts 
in lattice gauge theory is assumed throughout. It goes without saying that many 
interesting theoretical developments in lattice QCD are not covered and that the 
selected topics merely reflect the knowledge and preferences of the lecturer. 

Acknowledgements 

I would like to thank the organizers of the school, particularly Rajan Gupta, for the 
opportunity to give these lectures and for some very pleasant and stimulating weeks 
at Les Houches. Much of the material presented here is based on work done by the 
ALPHA collaboration. I wish to thank the members of the collaboration for their 
long-term commitment and an enjoyable team-work. 



2 Lattice effects and the continuum limit 

Nearly 20 years ago Symanzik set out to study the nature of the continuum limit in 
perturbation theory He showed that lattice theories can be described through 
an effective continuum theory in which the lattice spacing dependence is made ex- 
plicit. This work later led him to propose a method to accelerate the approach 
to the continuum limit, which is now known as the Symanzik improvement pro- 
gramme In this section we discuss the effective continuum theory for lattice 
QCD. A particular form of improvement will be considered in section 3. 

2.1 Preliminaries 

In the standard formulation of lattice QCD (which goes back to Wilson's famous 
paper of 1974 (0) the quark fields ip(x) reside on the sites x of the lattice and carry 
colour, flavour and Dirac indices as in the continuum theory. The gauge field is 
represented through a field of SU(3) matrices U(x, /x) where fi = 0, . . . , 3 labels the 
space-time directions. None of these is special since we are in euclidean space, but 
fi = is conventionally associated with the time direction. 

Under a gauge transformation A(x) the fields transform according to 

ip(x) -» A(x)ip(x), U(x, fj) -> A(x)U(x, (j,)A(x + a/i) -1 , 
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where fi denotes the unit vector in direction /i. The forward difference operator 

VMx) = -[U(x,p)ip(x + afi) - ip(x)} 
a 

is hence gauge covariant. In the same way a backward difference operator V^j may 
be defined and the lattice Dirac operator is then given by 

£ = H7M(V£ + V M )-aV*V M }. (1) 

The second term (which is referred to as the Wilson term) is included to avoid 
the well-known "doubler" problem. It assigns a mass proportional to 1/a to all 
modes with momenta of this order while the low-momentum modes are affected by 
corrections vanishing proportionally to a. In particular, only these modes survive in 
the continuum limit. 

The QCD action on the lattice may now be written in the form 

S = \ E tr i 1 " U (P)} + ffl4 E ^){D + moMs) ( 2 ) 
% P x 

with go being the bare gauge coupling and uiq the bare quark mass. For simplicity 
we shall assume that the quark mass is the same for all flavours. The sum in the 
gauge field part of the action runs over all oriented plaquettes p and U (p) denotes 
the product of the gauge field variables around p. 

Exercise: List all continuous and discrete symmetries of the lattice action and write down 
the corresponding transformations of the fields 

2.2 Lattice effects in perturbation theory 

In lattice QCD one is primarily interested in the non-perturbative aspects of the 
theory. Perturbation theory can, however, give important structural insights and it 
has proved useful to study the approach to the continuum limit in this framework. 
A remarkable result in this connection is that the existence of the limit has been 
rigorously established to all orders of the expansion [|5L|j . 

The Feynman rules on the lattice are derived straightforwardly from the lattice 
action by introducing the gluon field A fl (x) through 

U(x,n) = exp{g aA tl (x)}, A„(x) = A°(x)\ a , 

and expanding in powers of go . One also needs to fix the gauge but this will not 
be discussed here (details can be found in ref. |g for example). Compared to the 
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usual Feynman rules, an important difference is that the propagators and vertices 
are relatively complicated functions of the momenta and of the lattice spacing a. In 
particular, at tree-level all the lattice spacing dependence arises in this way. 
A simple example illustrating this is the quark-gluon vertex 

= -g \ a {in-Hp + p)v + °( a2 )}- ( 3 ) 

It is immediately clear from this expression that the leading lattice corrections to 
the continuum term can be quite large even if the quark momenta p and p' are well 
below the momentum cutoff ir /a. Moreover the corrections violate chiral symmetry, 
a fact that has long been a source of concern since many properties of low-energy 
QCD depend on this symmetry. 

Lattice Feynman diagrams with I loops and engineering dimension u> can be 
expanded in an asymptotic series of the form 

oo I 

a~"J2J2c kj a h [ln(a)P'. 

k=0j=0 

After renormalization the negative powers in the lattice spacing and the logarithmi- 
cally divergent terms cancel in the sum of all diagrams. The leading lattice correc- 
tions thus vanish proportionally to a (up to logarithms) at any order of perturbation 
theory. 

Exercise: Work out the exact expression for the quark-gluon vertex and verify eq. (||) 
2. 3 Non-perturbative test for lattice effects 

Sizeable lattice effects are also observed at the non-perturbative level when calcu- 
lating hadron masses, for example. An impressive demonstration of this is obtained 
as follows. Let us suppose that the quark mass is adjusted so that the mass of the 
pseudo-scalar mesons coincides with the physical kaon mass. This sets the quark 
mass to about half the strange quark mass and one thus expects that the mass of 
the lightest vector meson is given by 

m v ~ m K * = 892 MeV. 

Computations of the meson masses using numerical simulations however show that 
this is not the case at the accessible lattice spacings (see fig. 1). Instead one observes 
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Figure 1 : Calculated values of the vector meson mass (full circles) and linear extrapolation 
to the continuum limit (cross). Simulation data from Butler et al. (GF11 collab.) {jj. 

a strong dependence on the lattice spacing and it is only after extrapolating the data 
to a = that one ends up with a value close to expectations. 

The numerical simulations quoted here have been performed in "quenched QCD" 
where sea quark effects are neglected. It seems unlikely, however, that a very dif- 
ferent result would be obtained in the full theory. Perturbation theory and the 
calculations in quenched QCD rather suggest that large lattice effects already arise 
from the interaction between the gluon field and the valence quarks and these are 
not suppressed when the sea quarks are included. 

2.4 Effective continuum theory 

In phenomenology it is well-known that the effects of as yet undetected substructures 
or heavy particles may be described by adding higher-dimensional interaction terms 
to the Standard Model lagrangian. From the point of view of an underlying more 
complete theory, the Standard Model together with the added terms is then a low- 
energy effective theory. 

A similar situation occurs in lattice QCD, where the momentum cutoff may be 
regarded (in a purely mathematical sense) as a scale of new physics. The associated 
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low-energy effective theory is a continuum theory with action |2],[|] 

S eff = J d 4 x {C (x) + aU{x) + a 2 C 2 (x) + ...}, (4) 

where jCq denotes the continuum QCD lagrangian and the C^s, k > 1, are linear 
combinations of local operators of dimension 4 + k. The dimension counting here 
includes the (non-negative) powers of the quark mass m by which some of the fields 
may be multiplied. From the list of all possible such fields only those need to be 
considered which are invariant under the symmetries of the lattice theory. Many 
terms can also be eliminated using partial integration. 

In the following we shall be mainly interested in the lattice corrections of order a. 
All other corrections are proportional to higher powers of a and are thus expected 
to be less important at small lattice spacings. Taking the symmetries of the lattice 
action into account, it is not difficult to show that the effective lagrangian L\ may 
be written as a linear combination of the fields 

Ox = tfr icr^ v F^ u 7p, 

e> 3 = mtr {F^F^} , 

04 = m {ipry^D^ip - ipD^tp} , 

05 = m 2 tpip, 

where F^ u denotes the gauge field strength and the gauge covariant partial 
derivative. 

In the effective continuum theory, renormalized lattice fields are represented 
through effective fields of the form 

<PcS = 4>o + a<t>i + a 2 <p2 + ■■■ (5) 

The fields (fro, (fri, ■ ■ ■ which can occur here must have the appropriate dimension and 
transform under the lattice symmetries in the same way as the lattice field which is 
being represented. As an example let us consider the isovector axial current 

A»(x)=^) 7m75 |t>(x). (6) 

In this case the leading correction (fri can be shown to be a linear combination of 
the fields 

(0 6 )^ = r i\)^\T a a ilu D v i\) - ifr D u a^j 5 ^T a ifr, 
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(0 7 )« = ^{y> 7 5^>} 



For other low-dimensional fields, such as the vector currents and the scalar and 
pseudo-scalar densities, the general form of the 0(a) corrections is also easily found. 
The number of terms is however rapidly increasing when one considers operators 
involving four quark fields or higher derivatives. 

Exercise: Prove that any local field of dimension 5 respecting the symmetries of the lattice 
theory can be written as a linear combination of the basis 0\ , . . . , O5 up to derivative terms 

2. 5 Scope of the effective continuum theory 

So far we did not specify how precisely the lattice theory is related to the effective 
continuum theory. One must be a bit careful here, because additional lattice effects 
can arise when integrating correlation functions over short-distance regions. Such ef- 
fects are characteristic of the correlation functions considered and are not accounted 
for by the effective action and the effective fields. In the following attention will 
thus be restricted to position space correlation functions at non-zero (physical) dis- 
tances. This restriction is not very severe, since most physical quantities, including 
the hadron masses and matrix elements of local operators between particle states, 
can be extracted from such correlation functions. 

So let us consider some local gauge invariant field 4>{x) constructed from the 
quark and gauge fields on the lattice. For simplicity we assume that <j)(x) does not 
mix with other fields under renormalization. We then expect that the connected 
renormalized n-point correlation functions 



have a well-defined continuum limit, provided the renormalization factor is chosen 
appropriately and if all points kept at non-zero distances from one 

another. 

In the effective continuum theory, the lattice correlation functions are represented 
through an asymptotic expansion 



G n (x 1 ,...,x n ) = {Z$) n . . .4>{x n )) 



con 



G n (x!,...,X n ) = ((/>o(xi)...0 O (x n )) 



con 




con 
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n 

+(*Y1 (M x i) ■ ■ ■ M x k) ■ ■ ■ M x n)) con + 0(a 2 ), (7) 
k=l 

where the expectation values on the right-hand side are to be taken in the contin- 
uum theory with lagrangian Co. The second term is the contribution of the 0(a) 
correction in the effective action. Note that the integral over y in general diverges at 
the points y = Xk- A subtraction prescription must hence be supplied. The precise 
way in which this happens is unimportant, because the arbitrariness that one has 
amounts to a local operator insertion at these points, i.e. to a redefinition of the 
field 4>i(x). 

It should be emphasized at this point that not all the dependence on the lattice 
spacing comes from the explicit factors of a in eq. (0). The other source of a- 
dependence are the operators 4>i an d £k> which are linear combinations of some 
basis of fields. While the basis elements are independent of a, the coefficients need 
not be so, although they are expected to vary only slowly with a. In perturbation 
theory the coefficients are calculable polynomials in ln(a). 

2. 6 Elimination of redundant terms 

The list of operators given above for the 0(a) lattice corrections can be reduced using 
the field equations. Formal application of the field equations yields the relations 

Ox - 2 + 20 5 = 0, 4 + 20 5 = 0, (8) 

which allow one to eliminate 02 and 04, for example. 0% is also redundant since it 
can be expressed as a linear combination of 07 and 0g. 

There are two technical remarks which should be made here. The first is that 
in euclidean correlation functions the field equations are only valid up to contact 
terms. In eq. (|7|) these can arise when y gets close to one of the points x^. The 
possible contact terms there however amount to an operator insertion, which can be 
compensated by a redefinition of <j>i . The latter thus depends on exactly which basis 
of operators has been selected for the 0(a) effective action. Note that when the field 
equations are used to simplify <f>i, no contact terms can arise since all points Xk keep 
away from each other. 

The second observation is that renormalization and operator mixing must be 
taken into account when deriving operator relations from the field equations. In 
particular, the naive relations, eq. (||), are only valid at tree- level of perturbation 
theory. At higher orders they are replaced by linear combinations of all basis ele- 
ments with coefficients that depend on the coupling and the chosen normalization 
conditions. Barring singular events, it is then still possible to eliminate 02, 04 and 

06- 
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Exercise: Work out the free quark propagator on the lattice and check that its dependence 
on the lattice spacing matches with what is expected from the effective theory 



2.7 Concluding remarks 

Through the effective continuum theory a better understanding of the approach to 
the continuum limit is achieved. In particular, omitting terms that are redundant or 
which amount to renormalizations of the coupling and the quark mass, the leading 
lattice correction in the effective action is proportional to the Pauli term ip ia^F^tp. 
The lattice thus assigns an anomalous colour-magnetic moment of order a to the 
quarks. Very many more terms contribute to £2 and a simple physical interpretation 
is not easily given. The pattern of the lattice effects of order a? should hence be 
expected to be rather complicated. 



3 O(a) improvement 

While the general form of the effective continuum theory is determined by the field 
content and the symmetries of the lattice theory, the coefficients multiplying the 
basis fields in the effective action and the effective fields depend on how precisely 
the latter has been set up. With an improved lattice action and improved expressions 
for the local composite fields we may hence be able to reduce the size of the lattice 
effects and thus to accelerate the convergence to the continuum limit J^S] . The idea 
is similar to using higher-order discretization schemes for the numerical solution of 
partial differential equations, but the situation is not quite as simple because of the 
non-linearities and the associated renormalization effects. 

Improvement can be implemented in different ways and there is currently no 
single preferred way to proceed. The subject has been reviewed by Niedermayer [|j 
and many interesting contributions have been made since then. 0(a) improvement 
is a relatively modest approach, where one attempts to cancel all lattice effects of 
order a. This development started a long time ago [|-13] an d the method is now 



becoming increasingly usable although there are still a few loose ends. 
3.1 0(a) improved action 

In the following our aim is to construct an improved lattice action by adding a 
suitable counterterm to the Wilson action eq. (||). The counterterm should be chosen 
such that the order a term in the action of the effective continuum theory is cancelled. 
We are only interested in position space correlation functions of the type discussed 
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Figure 2: Graphical representation of the products of gauge field variables contributing to 
the lattice field tensor [eq. ([H])]. The point x is at the centre of the diagram where all loops 
start and end. 



above and may hence assume that the effective lagrangian C\ is a linear combination 
of the reduced basis of fields Oi, O3 and O5. It is then quite obvious that L\ can 
be made to vanish by adding a counterterm of the form 

8S = a 5 ]T [cxOtix) + c 3 6 3 (x) + c 5 5 (x)} , 

X 

where Ok is some lattice representation of the field Ok . 

Apart from renormalizations of the bare parameters and adjustments of the 
coefficients Ck, the discretization ambiguities that one has here are of order a 2 . In 
particular, we may choose to represent the fields tr {F^F^} and iptp by the pla- 
quette field and the local scalar density that already appear in the Wilson action. 
The terms proportional O3 and O5 then amount to a rescaling of the bare coupling 
and mass by factors of the form \+0{am). We ignore this rescaling for the time being 
and return to the issue later when we discuss the renormalization of the improved 
theory. 

For the improved action we thus obtain 

S = <Swii son + a 5 Y^ c sw ip(x)ja flu F flu (x)ip(x), (9) 

X 

where F^ v is a lattice representation of the gluon field tensor. Recall that the 
product of the gauge field variables around a plaquette in the (/x, i^)-plane is equal 
to 1 + a?F^ v + . . . in the classical continuum limit. A symmetric definition of the 
lattice field tensor is hence given by 

^0*0 = ^2 iQA x ) - Q»n(x)} (10) 

with Qfj, u (x) being the sum of the plaquette loops shown in fig. 2. 
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The improved action, eq. (pi), first appeared in a paper by Sheikholeslami and 



Wohlert in 1985 [12|. It did not receive too much attention at the time, because 
systematic studies of lattice effects were not feasible with the available computers. 
The situation has now changed and there is general agreement that improvement is 
useful or even necessary, particularly in full QCD where simulations are exceedingly 
expensive in terms of computer time. 

To achieve the desired improvement, the coefficient c sw multiplying the 0(a) 
counterterm in the improved action should be chosen appropriately. For various 
technical reasons this is a non-trivial task and the tuning of c sw has long been 
considered a weak point of the method. Sheikholeslami and Wohlert calculated the 
coefficient at tree-level of perturbation theory and Wohlert (l]| later obtained the 
one-loop formula 

c sw = 1 + 0.2659 xg 2 + O(g%). (11) 
More recently a strategy has been developed which allows one to compute c sw non- 



perturbatively using numerical simulations [14|. Chiral symmetry plays an impor- 
tant role in this context and the further discussion of this issue is hence postponed 
to section 5. 

Exercise: Write down the expression for the lattice field tensor explicitly and check that 
it transforms under the lattice symmetries in the expected way 



3. 2 Improvement of local fields 

Through the improvement of the action one achieves that on-shell quantities such 
as particle masses and scattering matrix elements approach the continuum limit 
with a rate proportional to a 2 . As can be seen from eq. @ , some terms of order a 
however remain uncancelled in correlation functions of local fields. For a complete 
improvement of such correlation functions one also needs to improve the fields and 
to specify appropriate renormalization conditions. 

Improved fields are constructed by adding local counterterms to the unimproved 
fields in such a way that the order o term in the associated effective fields is cancelled. 
The counterterms are linear combinations of a basis of operators with the appropriate 
dimension and symmetry properties. In the case of the isovector axial current, 
for example, the basis consists of a lattice representation of the fields Oq, Oh and 
Og. The first of these is redundant and may be dropped without loss. The term 
associated with Os amounts to a renormalization of the current and is hence only 
meaningful when a renormalization condition has been imposed. Since we have not 
done so at this point, we drop this term too and shall address the issue again in the 
next subsection. 
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We are then left with a single term and the improved axial current becomes 



(4)£ = ^ + ca4,P°, (12) 

where denotes the average of the forward and backward difference operators and 

P a (x) =^(x) 75 ±t<V(x) (13) 

the isovector pseudo-scalar density. The coefficient ca depends on oo and should 
be chosen so as to achieve the 0(a) improvement of the correlation functions of the 
axial current. In perturbation theory one obtains [15] 

c A = -0.00756 x gl + 0(g*) 

and in quenched QCD the coefficient is now also known non-perturbatively Eq|. 



Exercise: Show that the pseudo-scalar density, eq. (13), is already improved, i.e. no 0(a) 
counterterm is needed in this case 



3.3 Renormalization 

When deriving the improved action and the expression eq. ( |i~2|) for the improved 
axial current, all terms that amount to renormalizations of the coupling constant, 
the quark mass and the current by factors of the form l + 0(am) have been dropped. 
The reason for this has been that such factors should better be discussed in the 
context of renormalization which we now do. 

In QCD with only light quarks it is technically advantageous to employ so- 
called mass- independent renormalization schemes. If we ignore O(a) improvement 
for while, the renormalized coupling and quark mass in such schemes are related to 
the bare parameters through 

m R = m^Cflg.a/i), m q = m -m c . 

The important point here is that the renormalization factors only depend on the 
normalization mass \x and the coupling, but not on the quark mass. Note that an 
additive quark mass renormalization m c (equal to 1/a times some function of go) is 
required on the lattice, because chiral symmetry is not preserved. 

In the improved theory the definition of the renormalized parameters must be 
slightly modified, as otherwise one may end up with uncancelled cutoff effects of 



13 



order om q in some correlation functions. Prom our discussion above one infers that 
we should first rescale the bare parameters according to 



9o = 9q0- + b g am q ) , 
fh q = m q (1 + b m am q ) , 

and then define the renormalized coupling and mass through 

al = 9o Z g{9o,afJ>), 
m n = fh q Z m (gl,an). 

The coefficients 6 g and b m depend on g$ and must be chosen appropriately to cancel 
any remaining lattice effects of order am q . They have been worked out to one- loop 



order of perturbation theory [18,19| and b m is now also known non-perturbatively 
in quenched QCD |§. 

Factors of the form 1 + 0(am q ) must also be included in the definition of renor- 
malized improved fields. The renormalized axial current, for example, is given by 

(A R )«=^ A (l + 6 A am q )(A)£ 

and similar formulae apply to any multiplicatively renormalizable field. Even though 
the renormalized fields are written in this way, it should be noted that the b- 
coefficients do not depend on the renormalization conditions which are being im- 
posed. The latter should be set up at zero quark masses and only fix the Z-factors. 

Exercise: Consider the free quark theory and prove that b m = — \ in this case. Determine 
the renormalization factor needed to improve the position space correlation functions of the 
quark field 



3.4 Is 0(a) improvement effective? 

Numerical studies of lattice QCD using the improved action eq. (||) have been ini- 
tiated a few years ago [ 21-24] . The additional computational work required for 
the improvement term generally takes less than 20% of the total computer time 
spent for the simulations, also in the unquenched case [37|. An important question 
which one would like to answer with these studies is whether the lattice effects on 
the quantities of physical interest are indeed significantly reduced at the accessible 
lattice spacings. For the non-perturbatively determined values of c sw several col- 
laborations have recently set out to check this 1 25 — p0|] , but it is too early to draw 
definite conclusions. 
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For illustration let us again consider the calculation of the vector meson mass m v 
discussed in subsection 2.3. A preliminary analysis of simulation results from the 
UKQCD collaboration gives, for the 0(a) improved theory, m v = 924(17) MeV at 
a = 0.098 fm and m v = 932(26) MeV at a = 0.072 fm |7||||. These numbers do not 
show any significant dependence on the lattice spacing and they are also compatible 
with the value m v = 899(13) that one obtains through extrapolation to a = of 
the results from the unimproved theory (left-most point in fig. 1). Data at larger 
values of the lattice spacing recently published by Edwards et al. |3(| corroborate 
these findings. 

Other quantities that are being studied include the pseudo-scalar and vector 
meson decay constants and the renormalized quark masses. Compared to the hadron 
mass calculations, the situation in these cases is significantly more complicated, 
because one also needs to determine the appropriate renormalization factors and 
the improved expressions for the isovector axial and vector currents. The experience 
accumulated so far suggests that the residual lattice effects are indeed small if a < 0.1 
fm. Further confirmation is however still needed. 

3. 5 Synthesis 

At sufficiently small lattice spacings the effective continuum theory provides an 
elegant description of the approach to the continuum limit. Whether the currently 
accessible lattice spacings are in the range where the effective theory applies is not 
immediately clear, but the observed pattern of the lattice spacing dependence in 
the unimproved theory and the fact that O(a) improvement appears to work out 
strongly indicate this to be so. Very much smaller lattice spacings are then not 
required to reliably reach the continuum limit. It is evidently of great importance 
to put this conclusion on firmer grounds by continuing and extending the ongoing 
studies of O(a) improvement and other forms of improvement. 
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4 Chiral symmetry in lattice QCD 



Many properties of low-energy QCD can be explained through the flavour and chiral 
symmetries of the theory. For the analysis of the weak interactions of the quarks, 
these symmetries are also very important and the fact that chiral symmetry is not 
exactly preserved on the lattice thus is a source of concern. In this section we 
discuss how precisely the symmetry is violated and indicate how to make use of it 
even though it is only approximate. 

4-1 Chiral Ward identities 

We first consider the theory in the continuum limit and proceed formally, i.e. without 
paying attention to the proper definition of the correlation functions that occur. We 
assume that there is an isospin doublet of quarks with equal mass m and study the 
associated chiral transformations. 

In euclidean space the expectation value of any product O of local composite 
fields is given by the functional integral 

(0) = ±f Oe~ s . 

£ J fields 

Starting from this expression the general chiral Ward identity is obtained by applying 
an infinitesimal chiral rotation 

= U a ( X )±T a l5 lP(x), = UJ a {x)li,(x)^\T\ 

to the integration variables. In these equations r a denotes a Pauli matrix acting on 
the isospin indices of the quark field and u a (x) is a smooth function which vanishes 
outside some bounded region R. Since the Pauli matrices are traceless, it is obvious 
that the integration measure is invariant under the transformation and we thus 
conclude that 

(5SO) = (SO). (14) 

The chiral variation of the action and of the local fields of interest are easily worked 
out. In particular, for the variation of the action one finds 

5S = J d 4 xuj a {-dpA^ + 2mP fl } , (15) 

where and P a denote the isovector axial current and density [eqs. (|6|),([l3|)]. The 
derivation of this result involves a partial integration and it has been important 
here that Lo a (x) smoothly goes to zero outside some bounded domain as otherwise 
a boundary term might arise. 
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Figure 3: Choice of the region R when deriving the PCAC relation. The product of fields O 
is assumed to be localized in the shaded area away from R. 



By inserting particular products of fields O and by choosing different functions 
ui a (x), many useful relations (collectively referred to as the chiral Ward identities) 
are thus obtained. A few special cases will be considered in the following paragraphs. 

Exercise: Show that the axial current A^ and the isovector vector current Vjj; = ip^^\T a ij) 
transform into each other under infinitesimal chiral rotations 



4.2 PCAC relation 

Let us first assume that the fields in the product O are localized outside the region 



R as shown in fig. 3. The variation 50 vanishes under these conditions and eq. (14) 
reduces to (5S O) = 0. Since this equation holds for all functions to a (x) with support 
in R, it follows that 

{dpA'fflO) = 2m(P a {x)0). (16) 

At zero quark mass this identity expresses the conservation of the axial current and 
it is hence referred to as the partially conserved axial current (PCAC) relation. 

It should be emphasized that eq. fll6| ) holds irrespectively of the product O of 
fields as long as these are localized in a region not containing x. Of course this is just 
a reflection of the fact that the PCAC relation becomes an operator identity when 
the theory is set up in Minkowski space. Different choices of O then correspond to 
considering different matrix elements of the operator relation. 

Exercise: Derive the general Ward identity associated with the SU(2) flavour group. What 
happens if the quarks have unequal masses? 



4-3 Integrated Ward identities 

Another class of identities is obtained when O factors into two products of fields, 
Ojnt and O ex t, localized inside and outside the region R respectively (see fig. 4). In 
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Figure 4: Assumed localization regions of the field products 0- m t and ex t when deriving 
the integrated Ward identity eq. (|17|). 

the limit where ui a {x) goes to a constant in R, eq. ( |l~4| ) then becomes 
/ da, jl (x)u a (A a Jx)O int O ext ) = 

- (JO int Oext) + 2m f d 4 x w a (P a (x)O int Oext> , (17) 
JR 

where the integration on the left-hand side runs over the boundary of R. So if we 
set m = and insert O m t = A b u {y), for example, a little algebra yields 

/ da^x) (A a (x)A b v (y)O ext ) = ie ahc (V u c (y)O cxt ) (18) 
JdR P 

with Vy(y) being the isovector vector current. 

Eq. ( p!8[ ) is reminiscent of the well-known current algebra identities in Minkowski 
space. Integrating the axial current over the boundary dR effectively amounts to 
taking the commutator of the axial charge with the fields localized in R. In parti- 
cular, the Ward identity eq. flilf ) simply states that the commutator of the axial 
charge with the axial current is equal to the vector current. Similar Ward identities 
may be derived for the isovector vector charge and the whole set of current algebra 
identities is thus recovered in euclidean space. 

4-4 Euclidean proof of the Goldstone theorem 

At vanishing quark mass and if chiral symmetry is spontaneously broken through a 
non-zero value of {ipip), the Goldstone theorem asserts that there exists an isovector 
of massless bosons (the pions) which couple to the pseudo-scalar density P a (x). In 
the following lines we give quick proof of this theorem in order to illustrate the power 
of the chiral Ward identities derived above. 



18 



Our starting point is the PCAC relation 

(d,A°(x)P a m = o 

which holds for all x 7^ 0. Using Lorentz invariance we then deduce that 



{A a (x)P a (0)) = k- 



x 



K x 2 y 

for some constant k. Now if we insert O m t = P a (0) and O ex t = 1 in the integrated 
Ward identity eq. (|l7|), and choose R to be a ball with radius r centred at the origin, 
the relation 

/ da^x) (A*(x)P a (0)) = -| m 

J \x\=r 

is obtained. The constant k appearing above may thus be calculated and one ends 
up with 

3 — , v x^i 



(A«(x)P«(0)) = - — W,)-^ 



We can now see that the current-density correlation function is long-ranged if (ipip) 
does not vanish. In particular, the energy spectrum does not have a gap and the 
correlation function in fact has a particle pole at zero momentum. 

4-5 PCAC relation on the lattice 

We now pass to the lattice theory and first discuss the PCAC relation. One may be 
tempted to derive this identity following the steps taken in the continuum theory, 
but this would not work out because the chiral variation of the lattice action cannot 



be expressed in terms of the axial current and density alone [31]. 

We may, however, argue that the position space correlation functions of properly 
renormalized lattice fields should converge to the corresponding correlation functions 
in the continuum theory. In particular, if (A^)^ and (P R ) a denote the renormalized 
isovector axial current and density, we expect that 

(MA R )«(x) O) = 2m R ((P R ,r(x) O) + 0(a) (19) 

for any product O of fields located at non-zero distances from x (as before the symbol 
dfj, stands for the average of the forward and backward difference operators). 

The size of the error term in eq. ( |i~9| ) directly tells us how strongly chiral symme- 
try is violated on the lattice. In particular, in the 0(a) improved theory the error 
should be reduced to order a 2 if the coefficients c sw and ca have their proper values. 
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x -y +t 



Figure 5: Assumed shape of the region R in eq. j20|). The vector n M denotes the outward 
normal to the boundary dR. 



This may conversely be taken as a condition to fix these coefficients, a possibility 
which will be examined in more detail in section 5. 

Other applications of the PCAC relation in lattice QCD include the computation 



of the additive quark mass renormalization m c [16| and of the renormalized quark 



masses (see ref. [45] for a review and an up-to-date list of references). 



Exercise: Since the flavour symmetry is unbroken on the lattice, there exists an exactly 
conserved isovector vector current. Deduce the explicit expression for this current 



4-6 Chiral symmetry & renormalization 

A well-known theorem states that canonically normalized conserved currents need 
not be renormalized, because the possible values of the associated charge are fixed 
algebraically. On the lattice chiral symmetry is violated and the theorem conse- 
quently does not apply to the axial current. In general the renormalization constant 
Za is hence different from 1 and needs to be calculated. 

As has been noted a long time ago, Za may be determined through the chiral 
Ward identities [^|J3^]. A possible starting point is eq. (H) with R being the 
volume between two equal-time hyper-planes (see fig. 5). At zero quark mass and 
with properly normalized currents we then expect that 

a 3 ]T e abc n,((A R )«(x)(A R ) b u (y)O ext ) 
xedR 



2i <(y R )S(y)0«t> + 0(a) (20) 
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for any product O cxt of fields localized outside R. Note that all fields in these 
correlation functions are in space-time regions well separated from each other. In 
the 0(a) improved theory, the error term may hence be expected to be reduced to 
order a 2 . 

Since the left-hand side is proportional to the square of Z\, and since the axial 
current does not appear on the other side of eq. (^), it is quite clear that the 
renormalization constant can be determined in this way. We shall not go into any 
further details here, but the idea has been shown to work out in practice and in 
many cases Z\ is now known non-perturbatively |33|--j3rf| . 

Once the renormalization constant of the axial current has been determined, the 
chiral Ward identities may be used to relate the renormalization constants of fields 
that belong to the same chiral multiplet. The isovector pseudo-scalar density (Pr) 
and the isoscalar scalar density S R transform into each other under chiral rotations 
and are thus an example of such a multiplet. In this case we would replace A^(y) in 



eq. (20) by the scalar density. The chirally rotated field then appears on the right- 
hand side and by evaluating the correlation functions on both sides of the equation 
one is hence able to fix the relative normalizations of these fields. In principle this 
procedure generalizes to any chiral multiplet, but some care should be paid when 
operator mixing occurs. 

Exercise: On the lattice the local isovector vector current V® = il>j f j,^T a ip is not exactly 
conserved. How can one fix the normalization of this current? 



4-7 Concluding remarks 

Our discussion in this section shows that there are many possibilities to make good 
use of chiral symmetry in lattice QCD even though the symmetry is not preserved 
on the lattice. The assumption which has implicitly been made is that the lattice 
corrections to the chiral Ward identities are small at the accessible lattice spacings. 
In the unimproved theory this is not obviously the case, because chiral symmetry 
is violated at order a and these effects can be rather large as we have previously 
noted. Further impressive evidence for this will be given below and the bottom-line 
then is that 0(a) improvement (or maybe some other form of improvement) should 
better be applied if one is interested in calculating quantities that are sensitive to 
chiral symmetry. 
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5 Non-perturbative improvement 



As explained in section 3, the coefficients multiplying the various 0(a) counterterms 
in the improved theory must be tuned to some particular values to achieve the desired 
improvement. We now show that the coefficient c sw appearing in the improved 
action, eq. (|9|), can be determined non-perturbatively through the requirement that 
chiral symmetry is preserved up to terms of order a 2 . 

The correlation functions that we shall consider for this purpose are constructed 
in a finite space-time volume with Dirichlet boundary conditions in the time direc- 
tion. There are good reasons for this choice, but these will only become clear later. 
Many technical details will be omitted in the following to keep the basic argumention 
transparent. A complete treatment is given in refs. [14 lq|. 



5.1 Boundary conditions 

The lattices we are going to consider are assumed to have spatial size L and time-like 
extent T, with the time coordinate taking values in the range 0, a, 2a, . . . , T. We 
impose periodic boundary conditions in all space directions and Dirichlet boundary 
conditions at time xq = and xq = T. The dynamical degrees of freedom of the 
gauge field are thus the variables U(x, p) associated with the links in the interior of 
the lattice, while at the boundaries we require that 

U(x,k)\ XQ=0 = e aC , U(x,k)\ Xo=T = e aC \ 

for all k = 1,2,3. The matrix C is taken to be constant diagonal, 

. /0 a \ 
C = - 2 , 01 + 2 + 3 = o, 
\0 3 / 

and C is similarly given in terms of another set of angles 4>' Q . 

In the case of the quark field only half of the Dirac components are fixed at 
the boundaries. The quark propagator is then the solution of a well-posed bound- 
ary value problem. Explicitly, if we introduce the projectors P± = ^(1 ± 70), the 
boundary conditions are 

P + il)(x)\ X0=0 = p(x), P-i/)(x)\ X0 =t = p'(x), 

where p and p' are some externally given (anti-commuting) fields. The corresponding 
boundary conditions on the anti-quark field are 

ip(x)P-\ XQ=0 = p(x), ip(x)P + \ XQ = T = p'(x). 
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The field components i/j(x) and i(j(x) at times < xq < T remain unconstrained and 
represent the dynamical part of the quark and anti-quark fields. 

Exercise: Consider the free quark theory in the continuum limit with boundary conditions 
as specified above. Solve the Dirac equation for arbitrary smooth boundary values p, p' and 
prove that the solution is unique 



5.2 Correlation functions in finite volume 

The lattice action in finite volume is essentially given by the expressions that we 
have written down in the preceeding sections for the infinite volume theory. In the 
case of the 0(a) counterterm in eq. @, for example, the only modification is that 
the sum is restricted to all points x in the interior of the lattice. 

Expectation values of products O of local fields are now obtained in the usual 
way through the functional integral. The integration is performed at fixed boundary 
values, i.e. the integration variables are the unconstrained dynamical field variables 
residing in the interior of the lattice. In particular, the only dependence on the 
boundary values C, C , p, p, p' and p' arises from the action. 

An interesting option at this point is to include derivatives with respect to the 
quark and anti-quark boundary values in the field product O. The operator 

° a = - a *Y,T^M^-^- (21) 

is an example of such a product. In the functional integral the derivatives act on 
the weight factor e~ s and have the effect of inserting certain combinations of the 
dynamical fields localized near the boundaries of the lattice. More precisely, one 
finds that in nearly all respects the derivatives —5/5p(x) and 5/5p(x) behave like 
an anti-quark and a quark field at time xq = 0. 

For illustration let us consider the correlation function 

(A*(x)O a ) = {± [ A«(x)O a e- s ) . (22) 

I & "'fields J p =p= p '=p'=0 

Since one sums over the positions y,z in eq. (21), the operator creates a quark 



and an anti-quark at time zero with vanishing spatial momentum. The correlation 
function then is proportional to the probability amplitude that the quark anti-quark 
pair propagates to the interior of the space-time volume and that it annihilates at 
the point x (see fig. 6). As usual such pictures have an exact meaning if the quark 
lines are thought to represent quark propagators at the current gauge field fy 



23 



time 




space 



Figure 6: Sketch of the space-time manifold on which the lattice theory is set up. C and 
C are the boundary values of the gauge field. The irregular lines represent the trajectory 
of a quark anti-quark pair, which is created at time through the operator O a [eq. (|21|)]. 



5.3 How large are the chiral symmetry violations? 



In principle the error term on the right-hand side of the PC AC relation eq. (19) 
provides an estimate of the size of chiral symmetry violation in lattice QCD. The 
renormalization factors in the expressions for the renormalized improved axial cur- 
rent and the associated pseudo-scalar density, 

(A R )J = Z A (l + b A am q ){A« + c A a^P a }, 
(P R ) a = Zp(l + b P am q )P a , 

are however not known at this point and a straightforward calculation of the error 
term is hence not possible. 

Now let us define an unrenormalized current quark mass through 



m 



2(P a O a ) 

where O a is the operator introduced above. The PCAC relation then implies 

Zp(l + bpam q ) 



m 



Z A (1 + bAam q ) 



m R + 0(a). 



(23) 



(24) 
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The renormalization constants and the renormalized mass do not depend on the 
kinematical parameters such as the time x$ at which the axial current is inserted 
or the boundary values C and C of the gauge field. Changing these parameters 
basically means to probe the PCAC relation in different ways. So if we consider two 
configurations of the kinematical parameters at the same point (go, amo) in the bare 
parameter space, and if mi and 777-2 denote the associated values of m, it follows that 

777,1 — "72 = 0(a). (25) 

By calculating mi and m,2 we thus obtain a direct check on the size of the lattice 
effects in the PCAC relation. 

In the improved theory one expects that the differences in the calculated values 
of m are reduced to order a 2 . This may not be totally obvious, because improvement 
has only been discussed in the context of the infinite volume theory. The boundaries 
at xq = and xq = T indeed require the introduction of further 0(a) counterterms, 
but since the PCAC relation is locally derived, it can be shown that these affect the 



relation at order a only [14]. 



5.4 Tests of chiral symmetry at tree-level of perturbation theory 

At vanishing coupling and if C = C = 0, the quarks propagate freely and the 
current quark mass m is hence expected to coincide with the bare mass mo up to 
lattice effects. A straightforward calculation confirms this and one also finds that 
777 is independent of the time xq at which the axial current is inserted as it should 
be if lattice effects are negligible. 

For non-zero boundary values C and C the situation changes and large lattice 
effects are observed in the unimproved theory. To illustrate this we consider a 16 x 8 3 
lattice and set 

(01,02,03) = § (-7T,0,7r) , (01,02,03) = B (-57T,27T,37r) . (26) 

As shown in fig. 7 the corresponding values of m strongly deviate from the free 
quark value and they are also far from being independent of Xq. These effects 
almost completely disappear when the improved action is used with c sw = 1 (which 
is the proper value to this order of perturbation theory). Improvement thus works 
very well and it is possible to prove that the residual lattice effects seen in fig. 7 are 
of order a 2 as expected. 

The outcome of these calculations can be readily understood if we note that 
the boundary values C and C induce a background gauge field in the space-time 
volume, equal to the least action configuration interpolating between C and C. 
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Figure 7: Plot of am on a 16 x 8 3 lattice with amo = 0.01 at tree-level of perturbation 
theory. The boundary values of the gauge field are zero (open circles) or as given by eq. (p6|). 
Full circles and crosses correspond to c sw = and c sw = 1 respectively. 



For the boundary values chosen above, the background field can be shown to be a 
constant colour-electric field given by 

U(x,fx) = e aB ^ x \ B (x)=0, B k (x) = C + (C - C) x Q /T. (27) 

Now since the lattice effects of order a can be described by an effective Pauli term, 
their magnitude is proportional to the strength of the background field and the tests 
that we have performed are hence directly probing for these effects. 

Exercise: Check that the PCAC relation is violated by 0(a 2 ) effects only in the free quark 
theory on an infinite lattice 



5.5 Computation o/c aw 

Without improvement the lattice corrections to the PCAC relation are rather large 
also at the non-perturbative level. This is made evident through fig. 8, where the 
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Figure 8: Unrenormalized current quark mass m on a 16 x 8 3 lattice with spacing a = 0.05 
fm, calculated using the unimproved action. Open and full circles correspond to zero and 
non-zero boundary values of the gauge field. 



results of a numerical simulation of the unimproved theory in the quenched approx- 
imation are plotted. In this calculation the bare mass mo has been set to some 
value close to m c and the bare coupling has been chosen so that = 6.4, which 
corresponds to a lattice spacing of about 0.05 fm. Again one observes a strong 
dependence of the calculated values of m on the background field and also on xq. 

We may now attempt to adjust the coefficients c sw and ca so as to cancel these 
effects. This works remarkably well and it is possible to bring almost all points to 
a well-defined plateau (see fig. 9). Similar results are obtained for larger lattices, 
other boundary values and so on. In general improvement is very efficient, at least 
for a < 0.1 fm, and with properly chosen c sw and ca the remaining violations of 
chiral symmetry are small. 



5. 6 Results and further developments 

Proceeding in this way the coefficient c sw has been determined in quenched QCD 
|1£.17] and now also in full QCD with a doublet of dynamical quarks |3q] . In the 
case of full QCD the range of couplings covered extends from small couplings, where 
contact is made with the one- loop formula eq. (11), up to couplings corresponding to 
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Figure 9: Same as fig. 8, but using the improved action and the improved axial current with 
c sw = 1.60 and ca = —0.027. The horizontal line is drawn to guide the eye. 

6/^0 — 5.2, where the lattice spacing is around 0.1 fm according to some preliminary 
studies (see fig. 10). The non-perturbative values of c sw deviate considerably from 
the one-loop result at the larger couplings, which is the relevant range for calculations 
of the hadron masses and other other properties of these particles. 

The data shown in fig. 10 are well represented by the rational expression 



At the level of the numerical accuracy that one has been able to reach, this fit is 
no worse than the simulation data and it may at this point just as well be taken as 
the definition of c sw . An important remark in this connection is that different ways 
to calculate c sw , using various background fields etc., yield results that differ by 
terms of order a. There is nothing fundamentally wrong with this. The systematic 
uncertainty in c sw merely reflects the fact that O(a) improvement is an asymptotic 
concept valid up to higher-order terms. 

Once c sw is known non-perturbatively, one may attempt to calculate the coeffi- 
cients of other O(a) correction terms using similar techniques. So far these calcu- 
lations have been limited to quenched QCD and results have been published for ca 



-sw 



1 - 0.454 gl - 0.175 + 0.012 gg + 0.045 gg 
1 - 0.720 gl 



(28) 
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Figure 10: Simulation results for c sw in QCD with a doublet of dynamical quarks |3q|. The 
solid curve represents the fit eq. (E8h and the dotted line the one-loop formula eq. (|ll|) . 



16 1, the O(a) counterterm needed to improve the isovector vector current [39| and 
various ^-coefficients |36,20|. Some of these appear to be less accessible and new 
ideas will be required to be able to compute them pC|]. 



6 Non-perturbative renormalization 

One of the fundamental problems in QCD is to establish the connection between 
the low-energy sector and the perturbative regime of the theory. There are many 
instances where this is required, a particularly obvious case being hadronic matrix 
elements of operators, whose normalization is specified at high energies through the 
MS scheme of dimensional regularization. In the present section we address the issue 
in the context of lattice QCD and briefly describe some of the techniques that have 
been proposed to resolve it. 
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Table 1: Recent results for m s (quenched QCD, MS scheme at fi = 2 GeV) 



m s [MeV] 


reference 


122(20) 


Allton et al. (APE collab.) fjlj 


112(5) 


Gockeler et al. (QCDSF collab.) || 


111(4) 


Aoki et al. (CP-PACS collab.) Q 


95(16) 


Gough et al. Q 


88(10) 


Gupta & Bhattacharya Q 



1 Example 

Let us first discuss one of the standard ways to compute the running quark masses 
in lattice QCD. The need for non-perturbative renormalization will then become 
clear. Any details not connected with this particular aspect of the calculation are 
omitted. 

A possible starting point to obtain the sum m u + rn s of the up and the strange 
quark masses is the PCAC relation 

m l/ K = ( m u + m s ){0\ uj 5 s \K + ) 

(lattice effects are ignored here) . Since the kaon mass m K and the decay constant f K 
are known from experiment, it suffices to evaluate the matrix element on the right- 
hand side of this equation. On the lattice one first computes the matrix element of 
the bare operator (u7ss)i a t and then multiplies the result with the renormalization 
factor Zp relating (tt75s)i a t to the renormalized density u^y^s. 

Some recent results for the strange quark mass obtained in this way or in similar 
ways are listed in Table 1 (further results can be found in the review of Bhattacharya 
and Gupta ^5|). The sizeable differences among these numbers have many causes. 
An important uncertainty arises from the fact that, in one form or another, the 
one-loop formula 

Zp = 1 + £ {(2/?r) ln(a/x) +k}+ 0(5 ° ) 

has been used to compute the renormalization factor, where \x denotes the normal- 
ization mass in the MS scheme and k a calculable constant that depends on the 
details of the lattice regularization. Bare perturbation theory has long been known 
to be unreliable at the couplings of interest and without further insight it is clearly 
impossible to reliably assess the error on the so calculated values of Zp. Evidently 
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this difficulty would go away if we were able to compute the renormalization constant 
non-perturbatively. 



6.2 Hadronic renormalization schemes 

We now proceed to introduce the notion of a hadronic renormalization scheme and 
shall then be able to formulate the problem of non-perturbative renormalization in 
more physical terms. 

In a hadronic scheme the renormalization conditions are imposed at low energies 
by requiring a set of physical quantities, such as the hadron masses and suitable 
hadronic matrix elements, to assume prescribed values. We may, for example, take 
the pion decay constant f w = 132 MeV as the basic reference scale and fix the quark 
masses through the dimensionless ratios 



On the lattice the combination af w then becomes to a function of go alone and a 
variation in the coupling thus amounts to changing the lattice spacing at constant 
/tt (and fixed quark masses). 

The normalization of local fields can be fixed similarly by requiring the renor- 
malized fields O r to satisfy 



for some hadronic states | i } and | / }. If there is operator mixing further normaliza- 
tion conditions must be imposed to fix the finite parts of the mixing coefficients. 

After these preparations we can now clearly say what non-perturbative renormal- 
ization means. Basically what we require is to match a given hadronic scheme with 
the MS scheme of dimensional regularization (or any other perturbative scheme). 
This amounts to calculating the running coupling and quark masses at high energies 
as well as the renormalization constants Xq that one needs to convert from the 
renormalized operators C R to the corresponding fields the MS scheme. 

In principle such a computation is straightforward. The key observation is that 
all physical amplitudes are uniquely determined functions of the external momenta 
and the lattice spacing once a definite hadronic scheme has been adopted. We can 
then consider a set of correlation functions of local fields at large momenta and 
arrange the momentum cutoff ir/a to be far above these scales. Lattice effects are 
negligible under these conditions and by comparing the exact correlation functions 
with their perturbation expansion in the MS scheme one is hence able to extract the 
running parameters and renormalization constants at normalization scales /j, given 
in units of f n for example. 



If, 



(f\O n \i 



) = prescribed value 



31 



A subtle technical point here is that the MS scheme is only defined to every finite 
order of perturbation theory. The statement that the running coupling etc. can be 
calculated has, therefore, a precise meaning only in the limit fj,/f w — > oo where 
higher-order corrections become negligible. 

Exercise: How precisely are different hadronic schemes related to each other? Convince 
yourself that it suffices to solve the non-perturbative renormalization problem for one par- 
ticular hadronic scheme 



6.3 Meanfield improved perturbation theory 

One of the attempts to solve the non-perturbative renormalization problem goes back 



to an observation of Parisi [46], who noted that the often poor convergence of the bare 
perturbation expansion may be due to the appearance of so-called tadpole diagrams. 
Meanfield theory then led him to suggest that these diagrams can be resummed by 
replacing the coupling ao = (7g/47r through a P = chq/P, where P = 1 — ^g^ + . . . 
denotes the plaquette expectation value. Manipulations of this sort have later been 
discussed in greater detail by Lepage and Mackenzie [^?J and their recipes are now 
widely used. 

For illustration let us consider the series 

ce Ms(f 1 ) = a o + di(a/j,)al + d 2 (aiJ,)al + . . . , (29) 

which one obtains by calculating some correlation functions at large momenta using 
bare perturbation theory and comparing with the expansion in the MS scheme. For 
the standard Wilson action the expansion coefficients are known up to two loops 
|48-|50||, the result in quenched QCD being 



dl(a//) = 5.8836 — — ln(a/x), 

27T 

51 

d 2 (au) = [dUafi)} 2 + 8.7907 >y. 

4ir z 

Note that we do not require [i to be much smaller than the momentum cutoff ir/a. 
When deriving the expansion eq. (p9|), one only needs to know that the external 
momenta entering the lattice correlation functions that one has chosen to work out 
are well below this scale. The normalization mass, on the other hand, only appears 
in the continuum theory and can be arbitrarily large. 

To apply these formulae one selects some value of the bare coupling and de- 
termines the lattice spacing in the chosen hadronic scheme. This step involves the 
calculation of meson masses etc. and is hence performed using numerical simulations. 
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Knowing the lattice spacing in physical units and the bare coupling, the right-hand 
side of eq. (|29| ) can be evaluated straightforwardly for any given \x. On the presently 
accessible lattices, go is around 1 and 1/a is a few GeV at most. So if we set /U to 
some reasonable value, say 10 GeV, the one- and two-loop terms in the expansion 
make a contribution roughly equal to 35% and 17% of the leading term. These cor- 
rections are not small and there is, therefore, little reason to trust the results that 
one obtains. 

Now if we take Parisi's coupling a P as the expansion parameter, the terms in 
the resulting series tend to be significantly reduced. We may, for example, adjust 
the scale factor a^i so that the one-loop term vanishes and the two-loop correction 
in the resulting expansion 

a Ms(p) = a P + 2.185 x ctp + ... , ii = 2.633 x 1/a, 

then is indeed rather small (a P is around 0.13 at the bare couplings of interest). As 
a consequence one may be more confident in applying this form of the perturbation 
expansion, although it remains unclear how the truncation error can be estimated 
in a reliable manner. 



6.4 Intermediate renormalization 

An interesting method to compute renormalization factors that does not rely on 



bare perturbation theory has been proposed by Martinelli et al. [51|. The idea is 
to proceed in two steps, first matching the lattice with an intermediate momentum 
subtraction (MOM) scheme and then passing to the MS scheme. The details of the 
intermediate MOM scheme do not influence the final results and are of only practical 
importance. One usually chooses the Landau gauge and imposes normalization 
conditions on the propagators and the vertex functions at some momentum p. In 
the case of the pseudo-scalar density, for example, the renormalization constant 
Z P 10M is defined through 



7 MOM / 7 MOM 

Z 2 /Z'p x 



(30) 



where Z 2 MOM 



denotes the quark wave function renormalization constant and the dia- 
grams represent the full and the bare vertex function associated with this operator. 
On a given lattice and for a range of momenta, the quark propagator and the 



full vertex function can be computed using numerical simulations [51,52]. Z 2 
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and Zp OM are thus obtained non-perturbatively. The total renormalization factor 
relating the lattice normalizations with the MS scheme is then given by 

Zp(g , a/x) = Zp OM (g , ap)Xp OM (g Tm ,p/n), (31) 

with X^ OM being the finite renormalization constant required to match the MOM 
with the MS scheme. Xp OM is known to one-loop order of renormalized perturbation 
theory and could easily be worked out to two loops. 

While this method avoids the use of bare perturbation theory, it has its own 
problems, the most important being that the momentum p should be significantly 
smaller than n/a to suppress the lattice effects, but not too small as otherwise one 
may not be confident to apply renormalized perturbation theory to compute Xp OM . 
It is not totally obvious that these conditions can be met on the currently accessible 
lattices |53,|54 and careful studies should be made to check this. 



6.5 Recursive finite- size technique 

The computational strategies for non-perturbative renormalization discussed so far 
assume that all relevant physical scales can be accommodated on a single lattice 
which is sufficiently small for the required calculations to be performed using nu- 
merical simulations. As a consequence the energy range where the hadronic scheme 
can be matched with perturbation theory tends to be rather narrow and it is then 
not easy to control the systematic errors. 

As has been noted some time ago fl55[, this difficulty can be overcome by sim- 
ulating a sequence of lattices with decreasing lattice spacings. Any single lattice 
only covers a limited range of distances, but through the use of a finite-volume 
renormalization scheme it is possible to match subsequent lattices. Effectively one 
thus constructs a non-perturbative renormalization group. In a few steps it is then 
possible to reach much higher energies than would otherwise be possible. 

The technique is generally applicable and does not require any uncontrolled 
approximations to be made [|56|-59|. Perhaps the greatest difficulty is that the basic 



strategy does not easily disclose itself. One also needs to perform extensive numerical 
simulations to set up the non-perturbative renormalization group and there is usually 
a significant amount of analytical work to be done. 

In the following section we discuss the physical picture underlying the technique 
and introduce some of the key notions. The method itself will then be explained in 
section 8, taking the computation of the running coupling in quenched QCD as an 
example. 
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(a) (b) (c) 

Figure 11: Pictures illustrating various physical situations in finite volume, (a) Hadrons in a 
large volume, (b) a qq meson in a box of size L ~ 2 fm, and (c) quarks in the femto-universe. 



7 QCD in finite volume and the femto-universe 

In quantum field theory the physical information is encoded in the correlation func- 
tions of local operators and these are hence the primary quantities to consider. From 
statistical mechanics one knows, however, that certain properties of the system can 
often be determined more easily by studying its behaviour in finite volume. The 
calculation of critical exponents is a classical case where such finite-size techniques 
are being applied. 

The questions one would like to answer in QCD are not the same as in statistical 
mechanics, but the general idea to probe the system through a finite volume proves 
to be fruitful here too. In this section our aim mainly is to provide a qualitative 
understanding of what happens when the volume is decreased. Unless stated oth- 
erwise, periodic boundary conditions are assumed and the lattice spacing is taken 
to be much smaller than the relevant physical scales so that lattice effects can be 
ignored. 



7. 1 Physical situation from large to small volumes 

Let us first consider the case where the spatial extent L of the lattice is significantly 
greater than the typical size of the hadrons (box (a) in fig. 11). Single hadrons are 
practically unaffected by the finite volume under these conditions except that their 
momenta must be integer multiples of 2n/L. For multi-particle states the situation 
is a bit more complicated, because the particles cannot get very far away from 
each other. Two-particle energy eigenstates, for example, really describe stationary 
scattering processes. If there are no resonances the corresponding energy values 
differ from the spectrum calculated for non-interacting particles by small amounts 
proportional to 1/L 3 j6C 63] . 
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When L is 2 fm or so, a single hadron barely fits into the box and the virtual pion 
cloud around the particle may be slightly distorted (box (b) in fig. 11). In particular, 
via the periodic boundary conditions it is possible that a pion is exchanged "around 
the world" . One of the consequences of this effect is that the masses of the hadrons 
are shifted relative to their infinite volume values by terms of order e~ m ^ L [ B0| . 



At still smaller values of L, the quark wave functions of the enclosed hadrons 
are squeezed and one observes rapidly increasing volume effects. Eventually when 
L drops below the confinement radius, say L < 0.5 fm, the quarks are liberated and 
hadrons cease to exist (box (c) in fig. 11). QCD in such small volumes is referred 
to as the femto-universe |64]]. Evidently physics is very different in this little world 
from what one is used to. The chosen boundary conditions play an important 
role, for example, and the energy levels are separated by large gaps that increase 
proportionally to 1/L. 

Exercise: Work out the volume dependence of the particle mass in the massive A</> 4 theory 
to first order in A. Check that the volume effects vanish exponentially at large L 



7.2 Asymptotic freedom and the limit L — > 

Now if we go to the extreme and allow the box size to become arbitrarily small, we 
are entering the short distance regime of QCD and thus expect that the effective 
gauge coupling scales to zero according to the perturbative renormalization group. 

To make this a bit more explicit, let us consider the pure gauge theory on a lattice 
of infinite extent in the time direction and twisted periodic or Dirichlet boundary 
conditions in the space directions. The reason for this exotic choice is that perturba- 
tion theory is much more complicated if periodic boundary conditions are imposed 
due to the appearance of physical zero modes in the gluon action |65-|72^. In this 
theory an effective gauge coupling may be defined through 



--L/2 



where F q q(r, L) denotes the force between static quarks at distance r in finite volume. 
Note that L is the only external physical scale on which this coupling depends. 

Perturbation theory at finite L is more complicated than in infinite volume, but 
with the chosen boundary conditions no fundamental difficulties are encountered 
when deriving the Feynman rules. It is possible to set up perturbation theory 
directly in the continuum theory using dimensional regularization. Alternatively 
one may choose to work on the lattice and to convert to the MS scheme later (if so 
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desired). Proceeding either ways the expansion 

"qq = «ms(m) + [( n / 27r ) hx(pL) + k] a-^in) 2 + ... (32) 

may be derived and if we set fi = 1/L in this equation it is immediately clear that 
a q q behaves like a running coupling at scale L. 

In full QCD with any number quarks, essentially the same arguments apply and 
the important conclusion then is that QCD is perturbatively soluble at sufficiently 
small L. Eq. ( |32"| ) also provides a link between the femto- universe and the theory in 
infinite volume. In particular, if we manage to compute a qq - for some box size given 
in physical units, we could use the relation to determine the MS coupling at this 
scale. This may be surprising at first sight, but the equation only reflects the fact 
that the basic lagrangian is the same independently of the physical context. 

7.3 Finite-volume renormalization schemes 

It should be quite clear at this point that the box size L may be regarded as a 
reference scale similar to the normalization mass in the MS scheme. Renormalization 
schemes that scale with L are in fact easily constructed by choosing some particular 
boundary conditions and by scaling all dimensionful external parameters entering 
the renormalization conditions with the appropriate power of L. In particular, if 
the time-like extent T of the lattice is finite, the ratio T/L should be set to some 
definite value. 

The coupling a qq introduced above complies with this general definition and 
there are many ways to set up finite-volume renormalization conditions for the local 
operators of interest. The normalization of the isovector pseudo-scalar density, for 
example, can be fixed through 

« 6 E ((Pnr^)(Pn) a (y))\ Xo - yo=L/2 = constant, 

where the constant should be adjusted so that P R has the standard normalization 
at tree-level of perturbation theory. If zero modes are excluded by the boundary 
conditions, one can then show that 

P R = {1 + [-(2/tt) ln(/iL) + k] assin) + . . .} P m {») 

(with another constant k). From the point of view of perturbation theory, finite- 
volume renormalization schemes are hence as good as any other scheme. 

Exercise: Invent a finite- volume renormalization condition for the quark masses 
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Figure 12: Strategy to compute the running coupling, taking low-energy data as input and 
using the non-perturbative renormalization group to scale up to high energies. 



8 Computation of the running coupling 

Non-perturbative renormalization is required in many places in lattice QCD, but 
the calculation of the running coupling is certainly one of the most important cases 
to consider. As already mentioned in section 6, finite-size techniques allow one to 
trace the evolution of the coupling from low to high energies and we would now like 
to describe this in some more detail. Although the method is generally applicable, 
attention will here be restricted to the pure SU(3) gauge theory. This automatically 
includes quenched QCD, since the renormalization of the coupling is independent of 
the valence quarks, but in full QCD a separate calculation will be required. 



8. 1 Strategy 

The computation of the running coupling discussed in this section follows the arrows 
in the diagram shown in fig. 12, starting at the lower-left corner. In this plot the 
energy is increasing from the bottom to the top while the entries in the left and 
right columns refer to infinite and finite volume quantities respectively. 

In the first step the chosen hadronic scheme is matched at low energies with a 
suitable finite-volume renormalization scheme. The running coupling in this scheme 
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Figure 13: Definition of Sommer's scale r$. The dimensionless function r 2 F q q(r) is mono- 
tonically rising and ro is equal to the distance r where it passes through 1.65. 



is then evolved non-perturbatively to higher energies using a recursive procedure. 
Eventually the perturbative regime is reached and the A parameter can be deter- 
mined with negligible systematic uncertainty. The conversion to the MS scheme is 



trivial at this point (top line of fig. 12). 



8.2 Hadronic scheme 

Since we are considering the pure gauge theory, there are no quark mass parameters 
to be fixed and a hadronic scheme is defined simply by specifying a low-energy 
reference scale. We might take the mass of the lightest glueball, for example, but 
in practice this would not be a good choice, because glueball masses are difficult to 
compute accurately through numerical simulations. 

The force i*qq(r) between static quarks at distance r is more accessible in this 
respect and may be used to define a reference scale ro through p3| 



^F q q(r ) = 1.65 (33) 
(see fig. 13). Comparing with the phenomenological charmonium potentials, the 



number on the right-hand side of eq. ( p3[ ) has been chosen so that ro — 0.5 fm. 
One should however keep in mind that the pure gauge theory is unphysical and an 
assignment of physical units to ro should not be given too much significance. 

Exercise: At distances r > 0.3 fm, the available simulation data for the heavy quark 
potential are well represented through V qq (r) = or — j^- + c. Calculate the string tension o 
assuming ro — 0.5 fm 
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8.3 Finite-volume scheme 



We now need to specify the boundary conditions in finite volume and to define a 
renormalized coupling that scales with L. For a number of technical reasons, the 
coupling a q q introduced previously is not employed, but in principle this would be an 
acceptable choice. The definition given below is based on the Schrodinger functional 
and requires some preparation. 

Let us consider a lattice with boundary conditions as specified in section 5 (and 
no quark fields). Although we have not done so, it is possible to allow for arbitrary 
space-dependent boundary values and the associated partition function 



Z[C',C]= [ e-" 

./fields 



is then referred to as the Schrodinger functional. Using the transfer matrix formal- 
ism, it is straightforward to show that 2[C, C] is equal to the quantum mechanical 
transition amplitude for going from the field configuration C at time xq = to the 
configuration C at time xq = T Jft3-S4|. 



At small couplings the functional integral is dominated by the minimal action 
configuration U{x,fj) = exp{aB^(x)} with the specified boundary values. In partic- 
ular, the perturbation expansion of the Schrodinger functional, 

T[B] = - ]nZ[C, C] = \t [B] + Ti[B] + g 2 T 2 [B} + . . . , (34) 

% 

is obtained by expanding about this field, which in many respects plays the role of a 



background field. In eq. (34) the leading term is equal to the classical action of the 
background field, while the higher-order terms require the calculation of Feynman 
diagrams with propagators and vertices depending on B. 

The renormalizability of the Schrodinger functional has been studied in pertur- 
bation theory and it turns out that all ultra-violet divergencies can be cancelled 
by including the usual counterterms plus a few boundary counterterms in the ac- 



tion 1 79]. In pure gauge theories boundary counterterms in fact do not occur, be- 
cause there are no operators with the appropriate dimension and symmetries. The 
Schrodinger functional is hence a renormalized quantity in this case. In particular, if 
we choose C and C to depend on some parameter rj, it is clear that a renormalized 
coupling may be defined through [83] 

I' = f^/f} • (35) 

Note that L is the only external scale in this formula as it should be in a finite- volume 
renormalization scheme. 
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The precise choice of the boundary values is largely arbitrary and was made 
essentially on the basis of practical considerations |3(]|8j|. As in section 5 we take 
C and C to be constant abelian and set 

(01,02,0a) = \ (-1,0,1)tt -\ (-2,1,1)7?, 



' ' ' - i(-3,l,2)7r + |(-2,l,l)r 7 . 



3 V "! V " 2 
mjV2>V3; — \ ( — 3, 1, 2) 7T + 5 

For small 77 one can then prove that the background field eq. (|27|) is the unique 
absolute minimum of the action [83]. 

The calculation of g 2 through numerical simulations does not present any par- 
ticular problem. The important point to note is that 

dF/di] = (dS/dri) 

is an expectation value of some combination of the gauge field variables close to the 
boundaries. Since gluon zero modes are excluded through the boundary conditions, 
the commonly used simulation algorithms remain effective even when L is very small 
in physical units. 



Exercise: Show that the background field eq. (|27]) solves the lattice field equations. Are 
there any other gauge inequivalent solutions with the same boundary values? 



Matching at low energies 



Returning to the diagram in fig. 12, we first need to match the hadronic scheme with 
the chosen finite-volume scheme. In the present case this amounts to calculating the 
coupling g 2 at some box size L given in units of vq. Evidently the matching should 
be done at some convenient point, where both tq and g 2 can be accurately computed 
using numerical simulations. Preliminary calculations show that values of g 2 around 
3.5 are in this range. For definiteness we take 

Luygx = box size where g 2 = 3.480 

as the matching point and our task then is to determine the ratio Lmax/ r o- 

To this end let us choose some value of the bare coupling go where one is 
able to perform numerical simulations of physically large lattices. Using standard 



techniques, r^/a may then be computed with negligible finite- volume errors [73- 
fr^l - At the same coupling one can also simulate lattices with smaller sizes, say 
L/a = 6, 8, . . . , 16, and calculate g 2 in each case. The renormalized coupling is 
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0.1 0.2 



Figure 14: Simulation results for £ max /ro and linear extrapolation to the continuum limit 
(left-most point). 

monotonically rising with L/a and L m ax/o may hence be determined through inter- 
polation to the point where g 2 = 3.480 (if it is contained in the range of lattice sizes 
considered). 

Proceeding in this way the ratio 

L max /r = (imax/a)(ro/a) _1 

has been calculated at several couplings. Different values of go correspond to different 
lattice spacings and any variation in the results that obtains thus signals the presence 
of lattice effects. As shown in fig. [l4] the available data decrease linearly with the 
lattice spacing and after extrapolation to the continuum limit one gets 

^maxAo = 0.680(26). (36) 

L max is thus approximately equal to 0.34 fm in physical units. Note that the ob- 
served linear dependence on the lattice spacing is theoretically expected, even though 
there are no quarks here, because of the presence of the space-time boundaries (see 
refs. |33],|l4[] for further details). 

8.5 Non-perturbative renormalization group 

We now proceed to scale the renormalized coupling g 2 to high energies (small box 
sizes), taking its value at L = L max as initial datum. This will be achieved through 
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Figure 15: Extrapolation of g 2 (2L) to the continuum limit at various values of g 2 {L). 



a stepwise procedure in which L is decreased successively by factors of 2. 

Let us first consider the continuum theory, where the scale evolution of the 
coupling is governed by the renormalization group equation 



L^ = -P(g) = b g 3 + b 1 g 5 + ..., 



(37) 



with bo = ll(47r)~ 2 and b\ = 102(47r)~ 4 being the usual coefficients. This equation 
determines the coupling at any scale if it is known at some point. In particular, 
there exists a well-defined function a(u) such that 

g\2L)=a(g 2 (L)). 

a(u) may be regarded as an integrated form of the j3- function and is referred to as 
the step scaling function. 

The crucial observation now is that ct(u) can be computed on the lattice through 
numerical simulations. One first chooses some convenient lattice size, say L/a = 8, 
and adjusts the bare coupling so that g 2 (L) = u. In practice the tuning which is 
required here does not cause any problems. We may then increase the lattice size by 
a factor 2 and calculate g 2 (2L) at the same value of the bare coupling. The result 
of this computation is equal to a{u) up to lattice effects of order a which can be 
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1 1.5 2 2.5 u 



Figure 16: Simulation results for the step scaling function [^6 59 1. The curve is a polynomial 
fit with two parameters as described in the text. 

extrapolated away by repeating the calculation for several values of L/a. As shown 



in fig. |15| the simulation data are in fact independent of the lattice spacing within 
errors and the continuum limit can thus be easily reached. 

Proceeding in this way the step scaling function has been determined for a large 
range of couplings (see fig. 16). At the smaller couplings the data compare well with 
the expansion 

a(u) = u + squ 2 + s\u 3 + S2U + • • • , (38) 
which one derives from eq. (|37l) , the first two coefficients being 

s = 2 ln(2)6 , si = 2 ln(2)6i + ag. 

We may in fact fit all the data by truncating the series at some low order and using 
the coefficients S2, S3, . . . as fit parameters. 

Once the step scaling function is known, the evolution of the coupling can be 
computed straightforwardly by solving the recursion 

u = 3.480, u k = a(uk+i), k = 0,1,2,... 

The initial coupling uo corresponds L = -L max and one thus obtains the values of 
g 2 at L = 2~ k L mSLX . If we demand that a(u) should only be used in the range of 
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couplings covered by the data, the recursion can be applied 8 times at most. After 
six steps, for example, one gets 

f = 1.053(12), L = 2- 6 L max , 

where the error has been calculated by propagating the statistical errors to the fit 
polynomial and solving the recursion using this function. Careful studies have been 
made to check that this gives correct error bounds and that the results do not depend 
on the number of fit parameters once a good fit quality is achieved. 

Exercise: Derive a recursion relation for the coefficients in eq. ( |38| ) assuming the coef- 
ficients bk of the /^-function are known 



8. 6 Computation of the A parameter 

As shown in fig. [I?] the data for the coupling a = g 2 /4ir obtained in this way cover a 
large range of energies from approximately 600 lYleV to 150 GeV (using L m ax — 

0.34 

fm to convert to physical units). The perturbative scaling of the coupling sets in 
rather early and for a < 0.08 the data lie on top of the two- and three-loop curves. 
This may be a bit surprising, but one should not conclude that the absence of large 
corrections to the perturbative evolution is a general feature of the theory. In other 
schemes the coupling behaves differently in general and there is usually no way to 
tell in advance at which energy the non-perturbative contributions become small. 

We may now attempt to determine the A parameter by scaling the coupling all 
the way up to infinite energies and taking the limit 

A = lim j\{g 2 ), X(u) = (6 «r 6l/26 °e- 1/26 ° u . (39) 

This is actually not the best way to proceed, because the statistical and systematic 
errors would not be easy to trace. A better starting point is the exact expression 



A = -|-A(g 2 )exp ji J" du 



1 1 6i 

+ 



b{u) bou 2 b^u 
where b(u) is related to the /^-function through 

b(f) = -ma) = b g 4 + hg 6 + ... 



(40) 



Eq. (|40D may be proved straightforwardly by noting that the right-hand side is 
independent of L and that it has the correct behaviour for L — > 0. 
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Figure 17: Comparison of simulation results |56|5^] for the running coupling (data points) 
with perturbation theory. The solid (dashed) curve is obtained by integrating eq. (571), 
starting at the right-most point and using the 3-loop (2-loop) /3-function. 



To compute the A parameter through this formula, one selects some value of L, 
say L = 2~ 6 L max , where the coupling is known accurately and where it is sufficiently 
small that the perturbation expansion for the 6-function can be trusted. The series 
has recently been worked out to three loop s j85| , |86| and using this result to calculate 
the integral on the right-hand side of eq. (|40| ) one obtains 

A = 0.211(16)/L max . (41) 

Since the scale evolution of the coupling is accurately reproduced by perturbation 
theory at high energies, the systematic error arising from the neglected four-loop 
(and higher-order) corrections to the 6-function are expected to be small. Changing 
L from 2" 6 L max to 2" 8 L max in fact does not have any appreciable effect on the 
calculated value of the A parameter. The influence of a hypothetical four-loop 
correction can also be estimated by treating this term as a perturbation in eq. (|40|) 
and assuming an approximately geometrical progression of the expansion coefficients 
bk- The error bounds that one gets in this way are far below the statistical error 
and one again concludes that the higher-order contributions can be neglected. 
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Exercise: Estimate the systematic uncertainty on the A parameter which would arise if 
the three-loop coefficient b 2 would not be known 



8.7 Conversion to the MS scheme 

To one-loop order of perturbation theory the relation between the finite- volume and 
the MS coupling is given by |56] 



« = OmsO) + [(U/27T) lnOL) " 1-2556] a^( M ) 2 + 



A short calculation, starting from the definition eq. (|39|), then leads to the (exact) 
conversion factor A^g/A = 2.049 and if this result is combined with eqs. (|36| ) and 
([4l|) one ends up with 

A M 3 = 0.636(54)/r . (42) 

Note that all reference to the intermediate finite-volume scheme has disappeared at 
this point. Eq. ( f42| ) simply expresses the A parameter in the MS scheme in terms 
of the low-energy scale ro and thus provides the solution of the non-perturbative 
renormalization problem for the running coupling. In physical units, using ro = 0.5 
fm, one gets A^§ = 251(21) MeV, but one should not forget that the conversion to 
physical units is ambiguous in the pure gauge theory. The solid result is eq. (p^). 

The running coupling a-^ (fi) at normalization masses given in units of ro may 
finally be calculated by solving the perturbative renormalization group equation in 
the MS scheme, taking eq. Q42"| ) as input. The numbers that one obtains tend to be 
significantly lower than the experimentally measured values of the strong coupling 
constant, but since we have not included the quark polarization effects there is no 
reason to be worried. Sea quarks affect the evolution of the coupling (it becomes 
flatter) and they also influence the low-energy reference scale in some way which is 
difficult to foresee. 



9 Conclusions and outlook 

One of the greatest advantages of the lattice approach to QCD is that the theory is 
mathematically well-defined from the beginning. The problem of non-perturbative 
renormalization, for example, can be given a precise meaning in this framework and 
its solution then reduces to a mainly computational task. Ultimately the goal is to 
obtain a clean test of QCD by calculating several quantities with small statistical 
and systematic errors and comparing the results with experiment. This has not been 
achieved so far, but the theoretical advances described in these lectures allow one to 
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be more confident about the basic strategies. In particular, as far as the continuum 
limit is concerned it now seems unlikely that lattice spacings very much smaller than 
0.1 fm will be needed for a reliable computation of the hadron masses and decay 
constants. 

At this point the main problem in lattice QCD is that sea quark effects are not 
easily included in the numerical simulations. Although the available algorithms have 
improved over the past few years, the amount of computer time required to simulate 
even a moderately large lattice is still enormous. As a result of the rapid progress 
in computer technology, simulations of full QCD are however becoming increasingly 
feasible. In the meantime the quenched approximation may be used as laboratory 
to test new ideas and to study more complicated physical problems (K — > irn decays 
for example). 
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